Untwisting of a Strained Cholesteric Elastomer by Disclination Loop Nucleation 
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The application of a sufficiently strong strain perpendicular to the pitch axis of a monodomain 
cholesteric elastomer unwinds the cholesteric helix. Previous theoretical analyses of this transition 
' ignored the effects of Frank elasticity which we include here. We find that the strain needed to un- 

wind the helix is reduced because of the Frank penalty and the cholesteric state becomes metastable 
above the transition. We consider in detail a previously proposed mechanism by which the topolog- 
ically stable helical texture is removed in the metastable state, namely by the nucleation of twist 
disclination loops in the plane perpendicular to the pitch axis. We present an approximate calcu- 
lation of the barrier energy for this nucleation process which neglects possible spatial variation of 
^5 ' the strain fields in the elastomer, as well as a more accurate calculation based on a finite element 

! ! \ modeling of the elastomer. 

a . 

PACS numbers: 61.30.-v, 61.30.Jf, 61.41. +o 
Ch ' I. INTRODUCTION 

o ■ 

o ■ ... . 

i— — T Cholesteric liquid crystals are composed of chiral molecules that favor a helical twist in their equilibrium state. 
Because the helical pitch is typically comparable to the wavelength of visible light, cholesterics are useful as optical 
devices, e.g., as displays or in mirrorless lasing 0- Cholesteric elastomers have recently attracted theoretical 
and experimental interest for their novel optical and mechanical properties. Cholesteric elastomers consist of either 

■ intrinsically chiral polymers crosslinked to form a gel [3j or nematic polymers crosslinked in the presence of a chiral 
I solvent which when removed leaves an imprinted helical texture in the gel . As in the case of conventional cholesteric 

liquids, the helical pitch in the elastomers is comparable to the wavelength of visible light; as a result cholesteric 
elastomers possess a photonic band gap, which may be mechanically tuned to useful advantage in waveguiding and 
lasing applications 

The behavior of a cholesteric elastomer under mechanical strain has been studied both experimentally p| and 
theoretically El Q. This problem bears some similarity to that of a cholesteric liquid in a magnetic or electric 
field H,|H, E3- If a magnetic field is applied to a cholesteric liquid perpendicular to the cholesteric pitch axis, then 
as the field is increased the pitch increases and twist walls begin to appear with a spacing equal to the pitch. The 

■ pitch (and thus the separation of the twist walls) diverges at a critical value of the magnetic field. However, in the 
case of a cholesteric elastomer while twist walls also appear if a strain is applied perpendicular to the pitch axis, 
the helical texture imprinted at the time of crosslinking in an elastomer 0, 1 1 1| prevents long wavelength distortions 
that would increase the pitch. Instead, the material shrinks along the pitch axis thereby reducing the pitch. If the 

J> , Frank energy associated with gradients of the director field is ignored, then the cholesteric twist is eliminated at a 
critical value of the strain by abrupt reversals of the twist walls which are of zero width at the critical strain [(J Q • 
It was argued in Refs. 0, Q that the Frank energy can be ignored because in a typical cholesteric elastomer the 
length scale at which the Frank energy is comparable to the elastic energy is very small compared to the pitch. Thus, 
except for values of the strain very close to the critical value needed to unwind the helix, it is reasonable to ignore 
Frank elasticity. However, if the shear modulus of the elastomer is reduced or if the distribution of nematic polymer 
chains is made more isotropic, then Frank elasticity can have significant effects. In any case, independent of material 
parameters the Frank energy will become infinite as the width of the twist walls approach zero at the transition to 
the untwisted state. Thus, to fully understand the nature of the transition from the helical to the untwisted (i.e., 
nematic) state, the effects of Frank elasticity must be considered. In this paper we study theoretically the effect of 
Frank elasticity on this transition by minimizing the energy of a cholesteric elastomer including the Frank energy. 
Inclusion of the Frank energy leads not surprisingly to twist walls of finite width at the transition (and thus a lower 
value of the critical strain). Furthermore, Frank elasticity leads to the metastability of the helical state above the 
transition, raising the question of how the topologically stable helical twist is removed. We explore this question 
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by considering the mechanism first proposed in Refs. |a, |7| , namely, the nucleation of twist disclination loops in the 
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planes of the twist walls. A similar mechanism was analyzed by Friedel and de Gennes ^} m the case of a cholesteric 
liquid in a magnetic field. If the initial field strength is large enough so that the material is in the nematic state 
(i.e., the helix is unwound), then as the field is reduced to a value less than critical field, the nematic state becomes 
metastable, and the equilibrium cholesteric phase is nucleated by the creation of a twist disclination loop. Friedel and 
de Gennes carried out an approximate analytic calculation of the critical radius and activation energy of such a loop. 
We consider a similar approximate calculation here. The main drawback to this approximate calculation is that we 
assume that the strain field in the elastomer is spatially uniform. We improve upon this calculation by carrying out 
a finite element method (FEM) analysis which allows us to minimize the elastomer energy with respect to variations 
in both the director field and the rubber elastic degrees of freedom. 

This paper is organized as follows. In the next section we review the theoretical analysis of the unwinding transition 
in the absence of Frank elasticity, first discussed in Refs. @,[z|]. ln Sec. IIIII wc consider the effect of Frank elasticity 
on the transition, followed in Sec. HVI bv an analysis of the nucleation of twist disclination loops which eliminate the 
cholesteric twist. Concluding remarks are offered in Sec. [V] 

II. UNTWISTING OF THE CHOLESTERIC HELIX IN THE ABSENCE OF FRANK ELASTICITY 

The microscopic statistical-mechanical thoery of nematic rubber elasticity is a generalization of the classicalrubber 
free energy density to the case of nematics (cholesterics are locally nematic) 0] : 

Aubber= ^Tr^-^-r 1 'A) , (1) 

where fi is the rubber shear modulus, A is the gradient of the strain field: 

Xij = Sij + VjUi, i,j = x,y,z (2) 

where u is the displacement vector. The shape tensor £ Q corresponding to the director field no before the mechanical 
deformation is applied is given by: 

lo = 'ji+(*||-*i>ono (3) 
The inverse shape tensor t -1 in the presence of the applied strain when the director field is given by n is given by: 

i . n i 



The deformation tensor A is subject to the incompressibility condition det(A) = 1. 

Throughout this paper we consider the response of a cholesteric elastomer to a uniform strain, A, applied perpen- 
dicular to the pitch axis of the cholesteric. Assuming that the pitch axis is along the z direction, the director field in 
the absence of strain is given by, 

n = (cos 0o, sin 0o,O) (5) 

where 0o = irz/p and p is the pitch. After the deformation is applied, even in the presence of nucleated disclination 
loops, wc continue to assume that the director lies in the plane perpendicular to the pitch axis and can be written as, 



n(r) = (cos</>(r),sin0(r),O) (6) 

and 

Figure ^ shows schematically how the elastomer responds to small applied strains. Note that the directors at 
0o = 0, 7r and 0o = tt/2 do not change under strain. This situation holds as long as the imprinting of the cholesteric 
state before crosslinking dominates over the effects of Frank elasticity [llj], which we assume will the case here. 

As discussed in Ref. [6(, elastic compatibility requires that a uniform strain applied along the x-direction to an 
elastomer that is uniform in the x-y plane leads to the following strain tensor: 

(7) 
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FIG. 1: Nail-head schematic of one pitch length of a cholesteric elastomer (a) with no applied strain, i.e., A = 1, (an undeformed 
helix), and (b) in the presence of a uniform strain, A > 1. 



Note that this form assumes the elastomer is uniform in the x-y plane and need not be true if there is nonuniformity 
such as will arise when we consider nucleating disclination loops in Sec. IIVI The determination of the critical value 
of strain, A c , required to unwind the helix when Frank elasticity effects are neglected was carried out in Refs. 0,0]- 
Using Eqs. JHJl-Q the rubber energy density Eq. JQ) can be written as: 

/rubber = ^{(A 2 + \ 2 yy + r -^-[X 2 {rc 2 sl) + \ 2 yy (rc 2 s 2 - s 2 c 2 ) - 2XX yy (r - l) So c oS c]}, (8) 

where co,so,c and s are shorthand for cos 0o , sin 0o ; cos and sin0 respectively. Minimizing Eq. © with respect to 
4>{z) yields the algebraic equation: 

t an 26 = 2AA TO sin20 o 

9 (r - 1)(A 2 + A 2 y ) cos 20o + (r + 1)(A 2 - A 2 J ' 1 > 

where the chain anisotropy r = £\\/£±. The transverse strain X yy in Eq. is determined by minimizing the rubber 
energy density integrated over one pitch length with respect to X yy : 



d 



p 

dz /rubber ) = 0. (10) 



OXyy \JQ 

Initially as A increases from unity, is approximately equal to 0o- Moreover, 0(0o = qo z = 0) = and 0(0o = tt/2) = 
7r/2 because of the anchoring of the director to the elastic matrix; however, the latter condition will not be satisfied 
once the strain exceeds the critical value for unwinding. For 6 = tt/2 Eq. © has solutions = and = 7r/2, 
corresponding to the untwisted and twisted states respectively. The denominator in Eq. @ vanishes when 

2 / rXL - A 2 \ 

co s2 * = -l + — (pfr^), (ID 

and Eq. JHJl is then satisfied by = 7r/4. As long as C = r ^ v — A 2 > the director angle is approximately equal 
to 0o, leaving the pitch unchanged aside from the affine contraction \/\ X xX yy . However, as £ — > + , the value of 
0o satisfying Eq. Ijllfl approaches 7r/2 — indicating a twist wall with infinite dd/dz at 0o = 7r/2, i.e., a twist wall of 
infinitesimal width. 

For £ < 0, never attains the value 7r/4; i.e., the twist wall is removed. Therefore, C = defines the critical strain, 
A c , beyond which the twisted phase is no longer stable; when £ < 0, Eq. JJjJ yields 0(0o = tt/2) = 0. The critical 
strain A c is given approximately by r 2 ! 1 . For r = 1.9, e.g., A c w 1.23. 
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III. UNTWISTING IN THE PRESENCE OF FRANK ELASTICITY 



The principal effect of including Frank elastic energy in determining the untwisting of a cholesteric elastomer is 
that the twist walls will now have a finite width, with a thickness approximately given by the length scale £ ~ 
(l/(r — l))y/K/Jx 00) which is a measure of the length scale at which the contributions of the Frank and rubber 
energies to the total energy are comparable. In addition, Frank elasticity decreases the value of the critical strain 
A c , since if we imagine that the cholesteric state was imprinted by a chiral solvent removed after crosslinking, then 
the Frank energy penalty favors untwisting the cholesteric. Furthermore, as we shall demonstrate explicitly, Frank 
elasticity leads to the metastability of the twisted phase for A > A c , and the untwisted phase for A < A c ; thus, in the 
present context "critical" strain will refer to the value of the strain where the free energy of the twisted and untwisted 
states are equal. 

The total elastomer energy now consists of the rubber energy density Eq. JHJ and the Frank energy density fx (in 
the one elastic constant approximation): 



f K = \k{V4>) 2 



(12) 



Minimizing the sum of these two energies with respect to 4>{z) leads to the equation: 



dz 2 



■ gi(4> ) sm2(/) - # 2 (0o) cos 20, 



(13) 



where, 



9i W>o) = t(AL+AL ) cos 20o + - 



f\2 _ \2 \ 

y^xx yy' 1 



(14) 



and 



32 (0o) = -j^xx\ y sin 2</>o- 



(15) 



The boundary conditions accompanying Eq. (| 1 3D are: <fi(4>o = 0) = while 4>((f>o = 7r/2) = 7r/2,0, in the twisted and 
untwisted phases respectively. As in the previous section, the total energy of one pitch length of the cholesteric must 
also be minimized with respect to X yy : 



d 



dz (/rubber + /if) = 



(16) 



We solved Eqs. and Q16[) simultaneously using the shooting method choosing K = 1CP 11 J/m, [i = 10 5 
J/m 3 , r — 1.9, and p = 300nm The results are shown in Fig. [3 where the energies (integrated over half a pitch 
length which is sufficient by symmetry) of the twisted and untwisted states as functions of the applied strain are 
shown. 




/■p/2 

FIG. 2: Free energy F = — / dz(f Iuhhcl + f K ) vs. A for K/u = 10' 

Jo 



m and r = 1.9. 
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The energies of the twisted and untwisted states are equal when A c « 1.16, compared to the critical value A c ~ 1.23 
found in Refs. when the Frank energy is neglected (see Sec. [HI. When the Frank energy is neglected not only is 

the critical value of strain greater but the helical state is unstable above this critical value. In the presence of Frank 
elasticity the helical state remains metastable above the critical value of 1.16. 

Figure shows the dependence of X yy on A for the twisted and untwisted phases. Figure 0] shows the director 



Ayy 




FIG. 3: Xyy(X) for the untwisted and twisted states 








FIG. 4: The director solution <f>(4>o) f° r A = 1.2 in the twisted and untwisted states. Note that at this strain value the twisted 
state is metastable. 

solution <j) for the twisted and untwisted states at a strain slightly greater than the critical strain value of 1.16. Noting 
that the phase difference corresponding to the length scale £, A<j>o = n^/p, is of order 0.1 with our choice of material 
parameters, we estimate the width of the twist wall to be approximately 4£. Given that the Frank energy causes the 
twist walls to have finite width at the transition to the untwisted nematic state, it is necessary to ask how the twist 
stored in the helix is eliminated above the critical strain. 



IV. TWIST LOOP NUCLEATION 

At strains exceeding the critical value of 1.16 found in the previous section, the twisted state becomes metastable 
relative to the equilibrium untwisted state. As in Refs. 0,0 we consider the possibility that the decay of the twisted 
state and ensuing growth of the equilibrium untwisted state occurs via the homogeneous nucleation of twist disclination 
loops in the planes intersecting the z axis at z = (2n + l)p/2, n = 0, ±1, ±2, . . .., since this is where the Frank energy 
density is largest. Fig. |S] illustrates how the appearance and growth of a disclination loop in one of these planes leads 
to the removal of the helical twist. 
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FIG. 5: Schematic of the removal of twist in a cholesteric elastomer subject to strain A along the x axis via nucleation of a 
disclination loop: a) undeformed helix at zero strain, b) deformed helix in a strained elastomer with pinning at <f>o = 7r/2, c) 
removal of twist via the nucleation of a disclination loop on the plane corresponding to <j>o — tt/2. Note that the director fields 
in (b) and (c) differ substantially only over a small distance along the pitch axis of order £, above and below the disclination 
loop. 



A. Estimate of energy barrier for nucleation of disclination loop 



We first present an order-of-magnitude estimate of the energy cost of nucleating a circular disclination loop of 
radius R in one pitch length of a cholesteric elastomer subject to an applied strain slightly above the critical value for 
the untwisting transition. We assume that the strain field throughout the elastomer is given by the solution for the 
metastable twisted state (the upper curve on the right-hand side of Fig- • I n Sec. HVBl we relax this assumption 
and carry out a finite element method (FEM) calculation of the nucleation barrier. 

We use the approach of Ref. [Toj where the nucleation energy for a disclination loop in a cholesteric liquid in a 
magnetic field was estimated. The nucleation energy of a disclination loop in one pitch length is the difference in the 
total energy (Frank plus rubber elasticity) of the states with and without the loop: 



AE = 



J d 3 x (V(j)d) 2 + J d 3 xf, 



rubber 



rubber I <Pt 



(17) 



where the integrals are over one pitch length, and 4>d and <p t are the director phases in the presence and absence of 
the disclination loop respectively. Recall that in the absence of the loop the elastomer is in the metastable twist state; 
hence the "i" subscript on <f>. 

We assume first that <pd differs substantially from the uniformly twisted phase <j) t only within a distance £ above 
and below the plane of the disclination loop, z — p/2, and within a radial distance p = R + £. Furthermore, the 
leading contributions to the integration of (Vc^) 2 will come from z — > p/2. Therefore, we write 



'u = <Pt 



7T t'-p/2) , 

-e « h , 

2 \z-p/2 



p-R 



(18) 



where, by symmetry, we only need to consider z > p/2. Note that the exponential is needed in Eq. (|18f) as a 
convergence factor. 

In the limit z — > p/2 + , h ( ^Z^j 2 ) is a solution to Laplace's equation (since Frank energy dominates over rubber 
energy). In particular, 



p-R\ 1 1 / p-R 

77: =t: arctan — 

z-p/2) 2 7T \z-p/2 



(19) 



7 



Taking the asymptotic limit R 3> £, using Eqs. (|18(l and (|19|l . and assuming is constant over a distance £ about 
z = p/2, Eq. (|17|l becomes 



A£ » y Kflln ^ - tt 2 #(z = p/2)^i? 2 + y d 3 x [/ mbbe 



M - /rubber (^t)] 



(20) 



where a is the core size of the loop. The difference of the rubber energies appearing in Eq. l|2Uf) . which is nonzero only 
within a volume 2ttR 2 ^ about z — p/2, is readily evaluated using Eq. (jSJ) within our assumption of ignoring differences 
in X yy between the phases, and recalling that the disclination loop is located in the plane specified by tfio = tt/2. We 
find: 



/rubbered = 0) — /rubber(</>i = 7r /2) = 



n V A WU 



A 2 )- 



Therefore, 



AE 



-KRXn (f /a) - R 2 tt^/kJ1 



2.75 




E c 



(21) 



(22) 



where we have also included a core energy E cm& = 2iraKR (where a is a numerical factor of order one) and evaluated 
cf)' t (z = p/2) using the data in Fig. 0] with the numerical result 2.75ir/p. 

The nucleation energy increases linearly with R for small R, reaching a maximum for R = R c , where 



R a = 



\K\n (£/a)+2airK 



(23) 



and the corresponding nucleation energy (i.e., the nucleation barrier height) is given by, 



AE C = tt^/KJ/l 



2.75 — * / (rX 2 — A 2 ) 



R,. 



(24) 



Assuming the material parameters given in Sec. IIIII as well as choosing a 
corresponding value of X yy given by the upper curve in Fig. [3J we find: 



and, 



R c « O.lp w 30 nm, 



AE C w 10 5 K w 10" 6 J, 



O.Olp, and setting A = 1.2 (with the 

(25) 

(26) 



using values for p and K as in the previous section. Note that the energy scale of AE C is of order Kp ~ 10 5 K with 
our choice of material parameters. Recall that our calculation of AE C neglects any possible spatial variation in the 
strain field Ay, assuming that the strain is given throughout the material by the mean-field solution of Sec. IIIII for 
the metastable twist phase. In the next section we carry out a more accurate calculation of the nucleation energy 
allowing for the proper minimization of the energy with respect to the elastic degrees of freedom. 



B. Finite element calculation of the nucleation energy 



We improve upon the estimate obtained for the nucleation energy in the previous section by using the finite 
element method for elastic solids to minimize the total elastomer energy as a function of the displacement field u 
and the director field <p, subject to appropriate boundary conditions. By minimizing the energy with respect to the 
displacement field rather than the strain field Xij , we automatically satisfy the conditions of elastic compatibility |l6| . 

We minimize the energy in a rectangular parallelepiped bounded in the z direction by two planes containing 
disclination loops. The x and y dimensions of the parallelepiped are 2L, such that L/R >> 1. With a strain 
A imposed perpendicular to the z axis, symmetry implies that only one quadrant of the parallelepiped need be 
considered, e.g., the region specified by x > 0, y > 0, as shown in Fig. EI The finite element mesh shown in Fig. 
was chosen to be coarse in the z-direction to allow for manageable computation time while allowing a reasonably fine 
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mesh in the x — y plane. The smallest mesh spacing in the z-direction was chosen near the two disclination loops 
which are located on the bottom and top faces. Specifically, we set this mesh spacing to be 0.02p, which is less than 
£ ~ 0.05p, in order to capture the expected rapid variation of cj) in these regions. The core size a of the disclination 
loops was chosen to be a = 0.2£. 

The boundary conditions on the displacement field u are as follows (refer to Fig.[fjJ): on face (1), x = L, u x = (A— 1)L; 
on face (2), y = 0, u y = 0; on face (3), x = 0, u x = 0; on the bottom and top faces, Ui{x,y,z =p/2) = Ui(x,y,3p/2), 
i — x,y, while u z (x, y, z — p/2) — u z (x, y, 3p/2) + constant. No constraints were specified on face (4) in order not to 
bias the solution in any particular way. The boundary conditions on the director angle 4> are: on the bottom face, 
inside the disclination loop (j) = tt and outside 4> = 7r/2; on the top face, inside the loop 4> = it and outside <j) — 3tt/2. 
Finally, the position of the origin is fixed to suppress translation of the elastomer as a whole. 

The incompressibility of the elastomer is imposed numerically by including a term fg in the total energy density 
proportional to a bulk modulus B : 

f B = l/2B(&ct±-l) 2 , (27) 

with B 3> (J.- 

We minimized the total energy using the material parameters: A = 1.2, K/fi = 1CP 16 m 2 , r = 1.9, and B/fi = 10 5 . 
Figure |7| shows the solution for the director field cf>. From Fig. we estimate that the twist wall is approximately 
~ 0.25p « 5£, in agreement with our earlier estimate from Fig. 0] which was based on minimizing the total energy for 
the same applied strain but in the absence of disclination loops. We have verified that the size of the twist wall is 
insensitive to the mesh size in the z-direction; using sixteen elements in the z-direction gave a similar estimate for the 
twist wall width. To calculate the loop nucleation energy AE from our FEM solution we once again used Eq. JTJJ, 
excluding the core region from the integration and inserting the approximate value for the core energy, E c — 2nKR, 
used in Sec. IIV^I The nucleation energy as a function of the loop radius is shown in Fig.|Hlfor two different values 
of the bulk modulus B. The order of magnitude of both R c and Ai? nuc agrees very well with our earlier estimates, 
Eq. lIl'Mli and (|24|l . which neglected spatial variations in the strain field and assumed that the strain throughout the 
elastomer was given by the solution for the metastable twisted state. This agreement is not surprising if one examines 
the spatial variation of the transverse strain in our FEM solution, e.g., as illustrated in Fig. Elfor R — 2p. Note from 
the legend in the figure that the spatial variation in the strain across the mesh is less than 1%. Similarly, in Table [I] 
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0.50 0.56 0.61 0.67 0.72 0.78 0.83 0.89 0.94 1.00 1.06 1.11 1.17 1.22 1.28 1.33 1.39 1.44 1.50 



FIG. 7: The director angle 4>/n for A = 1.2, L/p = 5, R/p = 0.6, and B//J, — 10 . The main portion of the figure shows the 
variation of <f> near the upper disclination. The inset shows the entire x > 0, y > quadrant of the strained elastomer. 



we display the values of X yy , the spatial average of X yy over the mesh, and the accompanying standard deviations 
for a range of values of R. The standard deviations provide a measure of the spatial variation of the strain over the 
mesh, and we note that they are as in Fig. [5] less than 1% of the mean for all values of R listed. 

We also note from Table [I] that X yy decreases monotonically with increasing R, i.e., an increase in the fraction of 
the elastomer volume corresponding to the untwisted state which has a smaller value of X yy (see Fig. - This trend 
is consistent with our results from Sec. 11111 specifically, Fig. [21 which shows that X yy is smaller in the untwisted state. 
Thus, starting with no disclination loop in the metastable twisted state above the critical strain, we would expect 
that the average value of X yy would decrease as the loop begins to grow and a progressively larger fraction of the 
elastomer is occupied by the equilibrium untwisted state. 

We note from Fig.[5]that even at the edges of the parallelepiped X yy is less than its value 0.8805 in the twisted state 
in the absence of a disclination loop (see Fig.|3J). Thus, the fact that X yy decreases with increasing R does not seem 
to be due solely to the contribution of the untwisted portion of the elastomer within the loop, but instead indicates 
a nearly spatially uniform X yy that decreases with R. 

We have checked the accuracy of our FEM solution by minimizing the energy in the metastable helical state in the 
absence of a disclination loop using the same meshes as in Table [I] and obtain excellent agreement with our results 
from Sec. IIIII The standard deviations in Xu in the absence of the disclination loop do not differ substantially from 
those in Table |U suggesting that spatial variations in Xu are likely a result of the large but finite value of B / fi rather 
than a real difference between the presence or absence of the loop. The standard deviations tend to increase with 
decreasing B/n — they are about a factor of two larger for B/[i — I0 4 than for B//J, = 10 5 . 

We have also measured the shear strains, Ay, i ^ j, and found them to be very small, of order 10 -2 — 10 -3 , and 
decreasing in value with increasing B / /i. There is no readily discernible difference between their values in the uniform 
twisted state and the state with nucleated disclination loops. Recall that in an incompressible elastomer assumed to 
be uniform in the x-y plane (i.e., in the absence of a disclination loop) elastic compatibility requires that the shear 
strains vanish [|| (see Eq. (JJJ. 
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FIG. 8: AE vs. R/p for A = 1.2 and L/p = 5. (□): B/fj, = 10 5 . (o): B/fi = 10 4 



TABLE I: Mean values (based on spatial averaging over the entire FEM mesh with A = 1.2, B/fi — 10 5 ) of the diagonal strain 
components Xa for several values of the disclination loop radius R. The standard deviations in the data reflect the spatial 
variation of the strain over the mesh. 
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V. CONCLUSIONS 



We have explored theoretically the effect of a mechanical strain applied to a cholesteric elastomer perpendicular to 
the pitch axis, focusing primarily on the transition from the cholesteric (twisted) phase to the nematic (untwisted) 
phase. We have extended the analysis of previous researchers [fj, |7j by including Frank elasticity. Because the 
Frank energy penalizes director deformations about uniform alignment, there is a reduction in the magnitude of the 
strain needed to unwind the helix. Additionally, the penalty for a nonuniform director field causes the twist walls 
in the strained cholesteric elastomer to be of finite width at the unwinding transition, unlike the case where Frank 
elasticity is neglected and the transition occurs when the twist walls have zero width. Frank elasticity also leads to 
the metastability of the twisted state above the transition, prompting the question of the nature of the mechanism 
that transforms the twisted state to the untwisted one. To address this question we have adopted the proposal put 
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FIG. 9: Xyy contours for A = 1.2, L/p = 5, R/p = 2, and B/p = 10 5 



forth in Refs. [(J |jfl that the transition is driven by the nucleation of twist disclination loops in the planes of the 
twist walls. We have explored the consequences of this idea in two ways. First, following the work of Friedel and de 
Gennes |10| who considered a similar question in the case of cholesteric liquids in the presence of electric or magnetic 
fields, we analytically, though approximately, evaluated the energy cost of the disclination loops, and thus determined 
the critical radius and nucleation energy barrier. The main drawback to this approximate calculation is that we 
assume that the elastic strain field in the elastomer is uniform which cannot be absolutely correct if a disclination 
loop is present. We addressed this issue by carrying out a finite element method evaluation of the disclination loop 
energy which allowed us to properly minimize the energy of the elastomer with respect to both director and elastic 
degrees of freedom. This calculation produces results for the nucleation barrier and critical loop radius in very good 
agreement with the approximate calculation and shows explicitly that the strain field is very nearly uniform in the 
elastomer. 
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